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A  population  genetics  study  of  Anopheles  darlingi  (Diptera: 
Culicidae)  from  Colombia  based  on  random  amplified  polymorphic 
DNA-polymerase  chain  reaction  and  amplified  fragment  lenght 

polymorphism  markers 

Ranulfo  Gonzalez /+,  Richard  Wilkerson*,  Marco  Fidel  Suarez,  Felipe  Garcia, 
Gerardo  Gallego’*,  Fleiber  Cardenas,  Carmen  Elisa  Posso,  Myriam  Cristina  Duque“ 

Facultad  de  Ciencias  y  Facullad  de  Salud,  Universidad  del  VaJle,  Cali,  Colombia  ''Department  of  Entomology,  Walter  Reed  Army 
Institute  of  Research,  Silver  Spring,  MD,  US  *elntemational  Center  ofTropical  Agriculture.  Cali,  Colombia 

The  genetic  variation  and  population  structure  of  three  populations  of  Anopheles  darlingi  from  Colombia 
were  studied  using  random  amplified  polymorphic  markers  (RAPDs)  and  amplified  fragment  length  polymor¬ 
phism  markers  (AFLPs).  Six  RAPD  primers  produced  46  polymorphic  fragments,  while  two  AFLP  primer  com¬ 
binations  produced  197  polymorphic  fragments  from  71  DNA  samples.  Both  of  the  evaluated  genetic  markers 
showed  the  presence  of  gene  flow,  suggesting  that  Colombian  An.  darlingi  populations  are  in  panmixia.  Average 
genetic  diversity,  estimated  from  observed  heterozygosity,  was  0.374  (RAPD)  and  0.309  (AFLP).  RAPD  and 
AFLP  markers  showed  little  evidence  of  geographic  separation  between  eastern  and  western  populations; 
however,  the  FST  values  showed  high  gene  flow  between  the  two  western  populations  (RAPD:  F^  =  0.029;  Ntn: 
8.5;  AFLP:  F fr  =  0.051;  Nm:  4.7).  According  to  molecular  variance  analysis  (AMOVA),  the  genetic  distance 
between  populations  was  significant  (RAPD:d>ST  =  0.084;  AFLP =  0.229,  P  <  0.001).  The  FST  distances 
and  AMOVAs  using  AFLP  loci  support  the  differentiation  of  the  Guyana  biogeographic  province  population 
from  those  of  the  Choco-Magdalena.  In  this  last  region,  Choco  and  Cordoba  populations  showed  the  highest 
genetic  flow. 

Key  words:  amplified  fragment  length  polymorphism  -  Anopheles  darlingi  -  genetic  diversity  -  malaria  vectors  -  random  amplified 

polymorphic  DNA  -  polymerase  chain  reaction  -  Colombia 


Anopheles  darlingi  is  recognized  as  the  most  anthro- 
pophillic  and  endophagous  species  of  Anopheles  in  the 
Americas  (Fleming  1986),  the  primary  malaria  vector 
in  the  Neolropics  (Conn  et  al.  2001),  and  is  considered 
the  most  important  malaria  vector  in  the  Amazon  Basin 
(Conn  et  al.  1999).  It  is  most  frequently  found  in  the  hot 
and  humid  lowlands,  with  a  discontinuous  distribution 
from  Southern  Mexico  to  Northern  Argentina,  but  not 
reported  in  Costa  Rica,  Nicaragua,  and  Panama  (Lin- 
thieum  1988,  Manguin  et  al.  1999).  In  Colombia,  it  is 
found  up  to  an  altitude  of  450  m  in  the  premontane  zone 
of  the  Andes  (Fleming  1986),  which  is  characterized  by 
three  large  mountain  ranges  and  variety  of  geographic 
features  (Fig.  1)  that  clearly  separate  the  eastern  from 
the  western  regions  and  possibly  limit  the  free  gene  flow 
between  these  regions. 

Based  on  the  original  description  of  An.  darlingi  by 
Root  (1926),  phenotypic  variation  and  differences  in 
hematophagous  behavior  have  contributed  to  generating 
a  certain  degree  of  confusion  regarding  its  specific  taxo- 
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nomic  status  (Gal vaoetal.  1937,  Gal  vao&  Barreto  1938, 
Galvao  1940,  Linthicum  1988,  Rubio-Paiis  1998).  Given 
the  foregoing  and  the  great  importance  of  this  species,  a 
scries  of  population  studies  were  conducted,  using  mor¬ 
phological  and  molecular  tools  (Manguin  1999).  ori¬ 
ented  toward  clarifying  both  taxonomic  status  and  ge¬ 
netic  structure  within  its  range  of  distribution. 

Lounibos  and  Conn  (2000)  recently  reviewed  the  use 
of  molecular  markers  in  the  study  of  the  ecology  and 
genetics  of  An.  darlingi.  Its  biting  behavior  was  sum¬ 
marized  by  Rosa-Frcitas  et  al.  (1992),  and  Zimmerman 
(1992)  demonstrated  great  variability  in  post-blood  feed¬ 
ing/resting  behavior.  Depending  on  the  region  of  study, 
one,  two,  and  even  three  peaks  of  major  biting  activity 
have  been  observed,  as  well  as  variation  in  the  extension 
or  time  interval  of  peak  hematophagous  activity.  Hudson 
(1984)  and  Klein  and  Lima  (1990)  suggest  that  this  in¬ 
dicates  the  existence  of  a  species  complex.  Neverthe¬ 
less,  from  morphological  studies  did  not  reveal  substan¬ 
tial  differences  among  populations  throughout  the  area 
of  distribution,  including  populations  in  isolated  zones 
(Linthicum  1988,  Rubio- Palis  1998,  Manguin  et  al. 
1999).  According  to  Linthicum  (1988)  the  diagnostic 
characters  are  very  constant,  even  for  the  Guatemalan 
and  Honduran  populations  and  inlerpopulational  varia¬ 
tion  is  almost  equal  to  intrapopulational  variation. 

Additionally,  several  studies  on  this  species  have  shown 
great  genetic  population  variability  (Kreutzer  et  al.  1972, 
Steiner  et  al.  1982,  Tadei  et  al.1982,  Rosa-Freitas  et  al. 
1992,  Freilas-Sibajev  et  al.  1995,  Mendcs  dos  Santos  et  al. 
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Fig.  1 :  geographic  location  of  three  Colombian  Anopheles  darlingi  popu¬ 
lations:  Bet6  (Choco),  Tierralta  (Cordoba),  and  Granada  (Meta). 


1996,  Rafael  &  Tadei  1998,  2000,  1999,  Manguin  cl  al. 
1 999,  Conn  cl  aJ.  1 999,  Malafronlc  et  al.  1 999,  Conn  2001 ). 

Mendes  dos  Santos  et  al.  (1996)  analyzed  1 1  enzymes 
in  four  Amazonian  populations  and  found  variation  with 
respect  to  the  number  of  alleles  and  polymorphism  in  all 
the  studied  populations.  Later,  Mendes  dos  Santos  et  al. 
(1999),  analyzed  19  isozyme  loci  in  four  populations  from 
the  Amazon  River  Basin  and  found  that  they  were  geneti¬ 
cally  similar,  but  that  the  one  from  the  central  Amazon 
region  had  high  polymorphism  in  comparison  to  the 
marginal  populations,  which  were  primarily  monomorphic. 

Using  a  combination  of  morphological  data,  random 
amplified  polymorphic  DNA-polymerase  chain  reaction 
(RAPD-PCR),  isozymes  and  internal  transcribed  spacer 
(1TS2)  sequences  in  samples  from  Central  and  South 
America,  Manguin  et  al.  (1999)  found  that  although  the  frag¬ 
ments  produced  by  RAPD-PCR  showed  evidence  of  geo¬ 
graphic  partition,  all  populations  were  separated  by  rela¬ 
tively  short  genetic  distances.  According  to  Conn  et  al. 
(1999),  the  heterogeneity  observed  in  this  species'  biting 
behavior  (Rosa-Freitas  et  al.  1992),  variation  in  size 
(Lounibos  cl  al.  1995,  Charlwood  1996),  and  geographic 
differences  in  blood  seeking  periodicity  can  be  partly  ex¬ 
plained  by  population  structure  and  gene-flow  patterns  that 
could  also  affect  the  species'  vector  capacity. 

Analyzing  the  ITS2  sequences  from  five  Brazilian 
states,  Malafronte  el  al.  (1999)  found  that  they  were  al¬ 
most  identical.  However,  the  population  from  the  South¬ 
east  (Dourado,  state  of  Sao  Paulo)  had  4-5%  divergence 


with  respect  to  the  others.  This  is  in  contrast  with  find¬ 
ings  using  an  mtDNA  marker  (Conn  et  al.  1999),  in  which 
nine  collection  sites  were  compared,  including  Dourado 
in  Brazil,  with  no  substantive  differences  detected  in 
haplotype  or  nucleotide  diversity.  In  Colombia,  accord¬ 
ing  to  Gonzalez  (2001),  An.  darlingi  showed  morpho¬ 
logical  variation  and  different  biting  behavior  between 
populations  from  the  Choco- Magdalena  biogeographic 
region  and  Guyana. 

In  this  study,  based  on  observed  polymorphism,  us¬ 
ing  two  dominant  markers,  RAPD-PCR  and  amplified 
fragment  length  polymorphism  (AFLP),  we  analyzed 
three  populations  separated  by  a  geographic  gradient  in 
order  to  determine  whether  genetic  structure  exists  in 
Colombian  An.  darlingi  populations. 

MATERIALS  AND  METHODS 

Collection,  handling  and  identification  of  An. 
darlingi  -  The  analysis  was  carried  out  using  samples  of 
An.  darlingi  collected  from  three  locations  in  Colom¬ 
bia:  Bete  (06°  00'  00"  N  and  76°  46'  60”  W)  in  the  mu¬ 
nicipality  of  Medio  Atrato,  Chocd  Province;  Granada 
(03°  32'  19"  N  and  73°  W)  in  the  municipality  of  Granada, 
Meta  Province;  and  Tai  in  the  municipality  of  Tierralta, 
Cdrdoba  Province  (Fig.  I).  Mosquito  samples  were  col¬ 
lected  using  human  bait  following  standard  WHO  ( 1975) 
recommendations.  They  were  identified  morphologically 
(Linthicum  1988)  and  preserved  in  1.5  ml-micro- 
centrifuge  lubes  al  -70°C. 

DNA  extraction  and  RAPD-PCR  and  AFLP  condi¬ 
tions  -  DNA  isolation  was  done  according  to  Coen  et  al. 
(1982),  modified  as  in  Romans  (Black  IV  &  DuTcau 
1 997).  The  samples  were  analyzed  using  RAPD-PCR  and 
AFLP  markers. 

RAPD-PCR  -  Specific  amplification  reaction  condi¬ 
tions  were  standardized  in  order  to  obtain  well-defined, 
consistent  and  reproducible  banding  patterns:  1  ng//d 
DNA,  IX  buffer,  2.5  mM  MgCI2,  0.2  mM  dNTPs,  0.2 
mM  oligonucleotides,  1U  Taq,  adjusted  to  25  //I  with 
H20.  The  RAPD  amplification  was  done  with  a  seven- 
step  profile  (initial  denaturation  cycle  at  94°C  for  3  min; 
35  cycles  at  94°C  for  30  s,  35°C  for  45  s  and  72°C  for  1 
min  30  s.  The  final  extension  was  done  al  72°C  for  5 
min).  Followed  by  a  soak  temperature  of  4°C.  The  am¬ 
plifications  were  done  in  a  MJ  Research  PTC- 100  plate 
thermocycler.  Amplification  products  were  visualized  on 
1 .5%  agarose  gels  according  to  the  method  of  Sambrook 
et  al.  (1989).  A  total  of  75  RAPD-PCR  10-mer  oligo¬ 
nucleotides  from  Operon  Technologies®  were  first 
tested  on  three  individuals  from  each  locality.  Negative 
controls  were  used  in  all  reactions  in  order  to  detect 
artifacts  and  to  verify  the  amplification  reliability.  Six 
RAPD-PCR  primers  (AOS,  A 1 3,  B05,  B12,  B14,  W09) 
were  chosen  for  testing  64  samples  (24  from  Choc6,  26 
from  Meta,  and  14  from  Cordoba).  The  location  selec¬ 
tion  criteria  was  based  on  separation  by  distance  and  the 
presence  of  the  Andean  mountain  range  as  a  geographic- 
barrier  (Fig.  1).  Samples  were  run  in  gels  at  300  V  in  a 
horizontal  electrophoresis  chamber  (Life  Technologies 
Inc.)  for  approximately  1  h.  lambda  DNA,  digested  with 
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Fig.  2:  majority  consensus  tree  UPGMA  in  three  Colombian  Anopheles  darlingi  populations  (Cordoba,  Chocd.  and  Meta)  analyzed  from  1 97  amplified 
fragment  length  polymorphisms  loci.  The  branch  number  indicates  Ihe  times  that  the  topology  was  explained. 


Pst  I,  was  used  as  the  molecular-size  standard  in  each  run. 
The  banding  patterns  were  visualized  with  clhidium  bromide 
(0.8  /<g/ml).  Gels  and  banding  patterns  were  recorded  with 
EagleEye  II™  (Stratagcne)  equipment. 

AFLP  -  An  analytical  system  II  kit  (GIBCO  BRL)  was 
used.  Sixteen  AFLP  primer  combinations  were  iniatilly 
screened  and  two  combinations  were  chosen  from  those: 
(E-AC/M-CAA  and  E-AG/M-CAG).  Preamplification 
was  done  wilh  a  20-cycle  PCR  (30  s  at  94°C,  60  s  at 
56°C,  and  I  min  at  72°C,  soak  at  4°C).  The  radioactive 
labeling  of  primer  EcoR  1  (selected)  was  done  with 
dATP|y32P|.  The  selective  AFLP  amplification  was  done 


with  an  I  l-step  profile  (initial  cycle  at  94°C  for  30  s, 
65°C  for  30  s,  and  72°C  for  1  min;  12  cycles  at  94°C  for 
30  s,  65°C  for  30  s  (lowering  0.7°C  per  cycle  each  time), 
and  72°C  for  I  min;  and  23  cycles  at  94°C  for  30  s,  56°C 
for  30  s,  and  72°C  for  1  min).  Both  the  preamplification 
and  selective  amplification  were  done  in  a  MJ  Research 
PTC- 1 00  thermocycler.  Seventy  samples  (28  from  Choco, 
27  from  Meta  and  15  from  Cordoba)  were  tested. 

Amplification  products  were  separated  on  6%  poly¬ 
acrylamide  vertical  gels  at  constant  current  (~  100  W 
and  approximately  1800  V).  The  gel  was  placed  on 
Whatman  3  MM,  covered  with  vinyl  film,  dried  and  ex- 
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TABLE  I 


Average  heterozygosity  expected  in  three  Colombian  populations  of  Anopheles  darlingi,  using  three  methods  with  randon  amplified 
polymorphic  DNA-polymcrase chain  reaction  (RAPD-PCR)  and  amplified  fragment  length  polymorphism  (AFLP)  loci 


Marker 

Subpopulalion 

Average 

heterozygosity 

expected 

Bayesian  method 

Uniform  q 
distribution 

Nonuniform  q 
distribution 

Meta 

0.2835 

0.3024 

0.2990 

AFLP 

Choco 

0.2474 

0.2741 

0.2683 

Cordoba 

0.2326 

0.2824 

0.2725 

Total 

0.3034 

0.3093 

0.3089 

Meta 

0.3579 

0.3679 

0.3688 

RAPD-PCR 

Choco 

0.3534 

0.3706 

0.3745 

Cordoba 

0.3567 

0.3833 

0.3862 

Total 

0.3742 

0.3785 

03804 

TABLE  II 

Estimate  of  F^and  Nm  among  Colombian  Anopheles  darlingi  populations  with  random  amplified  polymorphic  DNA-polymcrase 
chain  reaction  (RAPD-PCR)  and  amplified  fragments  length  polymorphism  (AFLP)  loci 


Marker 

Population  pairs 

Wright 

Theta 

Lynch  &  Milligan 

]'ST 

Nm 

fST 

Nm 

Fsr 

Nm 

RAPD-PCR 

Cordoba-Meta 

0.055 

4.30 

0.090 

2.50 

0.071 

3.30 

Choco-Meta 

0.049 

4.90 

0.074 

3.10 

0.065 

3.60 

Cdrdoba-Chocd 

0.029 

8.50 

0.023 

10.80 

0.008 

30.40 

Total 

0.060 

3.90 

0.068 

3.40 

0.057 

4.20 

AFLP 

Cdrdoba-Mcta 

0.106 

2.10 

0.224 

0.90 

0.220 

0.90 

Choco  -  Meta 

0.111 

2.00 

0.224 

0.90 

0.213 

0.90 

Cdrdoba  -  Choco 

0.051 

4.70 

0.103 

2.20 

0.0% 

2.40 

Total 

0.127 

1.70 

0.199 

1.00 

0.200 

1.00 

TABLE  III 

Results  of  AMOVA  in  three  Colombian  Anopheles  darlingi  populations  with  random  amplified  polymorphic  DNA-polymerase  chain 
reaction  (RAPD-PCR)  and  amplified  fragments  length  polymorphism  (AFLP)  loci 


Mean 

Variance 

Source  of 

Sum  of  squares 

squares 

components 

Distances 

Significance 

Marker 

variation 

(SS) 

(MS) 

(%)° 

(Eyr) 

(P) 

RAPD-PCR 

Between 

55.739 

111.478 

8.43 

0.084 

<0.001 

populations 

Within 

19.508 

1190.015 

9137 

populations 

Between 

populations 

289.256 

144.628 

22.93 

0.229 

<0.001 

AFLP 

Within 

populations 

1259.747 

18.526 

77.07 

a  :  by  variance  components. 


posed  to  x-ray  film,  developed  and  kept  for  further  analy¬ 
ses.  Data  from  the  RAPD  and  AFLP  banding  patterns  from 
each  sample  were  coded  in  a  binary  data  matrix. 

Analysis  of  the  RAPD  and  AFLP  patterns  -  A  dis¬ 
tance  matrix  (1-M)  was  generated  with  the  RAPDPLOT 
program  of  Black  IV  (1993,  1995).  This  matrix  was  es¬ 
timated  from  the  fraction  of  pairings  (A/),  using  the  for¬ 
mula  M  -  NAffNT,  where  NAB  is  the  total  number  of  pair¬ 
ings  between  individuals  A  and  B  (for  both  the  absent  or 
present  bands)  and  the  total  number  of  “loci”  (frag¬ 
ments)  in  the  study.  These  were  then  used  to  generate 
the  respective  sets  of  100  matrices  read  by  PHY  LIP  3.5C 


(Fclscnstein  1993).  Each  matrix  of  the  set  was  collapsed, 
using  "Neighbor  joining"  algorithm  with  the  UPGMA 
method  for  generating  dendrograms  from  which  the  de¬ 
finitive  tree  was  obtained  by  means  of  the  strict-con- 
sensus  and  majority-rule  options  Black  IV  (1993,  1995). 

In  order  to  analyze  the  consistency  with  which  the 
RAPD  and  AFLP  datasets  supported  the  phenelic  ratio 
between  the  possible  taxa,  a  bootstrap  analysis  of  100 
pseudoreplicates  was  run,  using  the  RAPDBOOT  1.0 
program  (Black  IV  1995). 

Population  genetics  analysis  -  Genetic  variation  and 
the  genetic  structure  of  the  three  An.  darlingi  populations 
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were  assessed  by  calculating  heterozygosity,  the  FiT sta¬ 
tistics,  migration/generation  rates  (Nm),  molecular  analy¬ 
sis  of  variance  (AMOVA)  and  the  genetic  distances. 

Population  diversity  and  genetic  differentiation  were 
analyzed  using  the  allelic  RAPD  and  AFLP  allele  fre¬ 
quencies.  Three  methods  were  used:  the  traditional  square 
root  of  q  (Apostol  et  al.  1996).  the  Lynch  and  Milligan 
(1994)  method,  and  a  Bayesian  method  proposed  by 
Zhivotovsky  (1999).  Because  of  the  difficulties  of  ob¬ 
taining  heterozygosity  from  dominant  markers  (Black 
IV  1993,  1995,  Black  IV  &  Munstermann  1996,  Taba- 
chnick  &  Black  IV  1996,  Yan  et  al.  1999),  we  calculated 
it  according  to  Black  IV  and  Munstermann  (1996)  and 
Apostol  et  al.  (1996),  who  take  into  account  the  absent 
bands  (q)  and  estimate  of  the  null  allele  frequency  q/i) 

at  locus  /  (/=  1  ,...,L)  in  populations  j(j  =  1 . r)  as:  qfl) 

-  V,\j  (/),  where  jc-(i)  is  the  null  recessive  homozygolcs 
frequency  in  population  j  at  locus  i. 

The  Fst  values  and  the  Nm  migration  rates  were  cal¬ 
culated  using  the  RAPDFST  program  (Black  IV  1995), 
assuming  RAPD  locus  and  AFLP  dominance  and  a  popu¬ 
lation  at  Hardy- Weinberg  equilibrium.  The  RAPDFST 
program  estimates  the  F^  from  the  formula  proposed 
by  Wright  (1951).  It  also  estimates  the  effective  migra¬ 
tion  rate  (Nm)  and  the  population  structure  (0^-). 

Genetic  variation  components  within  and  among 
populations  were  estimated  by  AMOVA  using  W1NA- 
MOVA  1.5  (Excoffier  et  al.  1992)  and  AMOVA-PREP 
1.01  (Miller  1998)  programs.  The  <I>ST  value,  which  is 
analogous  to  represents  the  population  structure. 
Significance  levels  and  the  F  values  w'ere  computed  1000 
times  by  nonparametric  exchange  procedures. 

In  addition,  the  genetic  distances  between  the  popu¬ 
lations  were  analyzed  with  the  RAPDD1ST  1.0  program 
(Black  IV  1997),  Based  on  each  RAPD  or  AFLP  data  set, 
a  distance  matrix  (1-S)  was  generated  (Nei’s  similitude, 

1 972);  S  =  2Nab/(N a  +  NB),  where  NAB  is  the  number  of 
fragments  that  individuals  A  and  B  share  in  common,  NA 
is  the  number  of  fragments  from  individual  A  and  NB,  the 
number  of  fragments  from  individual  B.  The  consistency 
with  which  the  data  set  supported  the  estimated  ratios 
between  the  populations  was  evaluated  by  means  of  a 
bootstrap  of  100  pscudoreplicalcs  as  described  earlier. 

RESULTS 

In  the  three  populations,  the  six  RAPD  oligos  (A05, 
A 13,  B05,  B12,  B14,  W09)  produced  45  polymorphic 
bands,  while  the  two  AFLP  primer  combinations  (E-AC/ 
M-CAA  and  E-AG/M-CAG)  produced  197. 

Analysis  of  the  RAPD  and  AFLP  patterns  -  The  clus¬ 
ter  analyses  with  the  RAPD  markers  using  the 
RAPDPLOT  (Black  IV  1993),  based  on  two  genetic  dis¬ 
tance  matrices  (1-M  and  1-S)  and  the  subsequent 
bootstrapping,  showed  different  individuals  clustering 
degrees  from  the  three  populations  analyzed  (Choco, 
Meta,  and  C6rdoba)  (cluster  not  presented).  The  low  lev¬ 
els  of  probability  obtained  in  these  clusterings  with  AFLPs 
do  not  indicate  a  significant  geographic  partition  (Fig.  2). 

Population  genetics  analyses  -  Both  the  RAPD  and 
AFLP  patterns  revealed  great  genetic  diversity  (Nei 
1972).  The  expected  average  heterozygosity  is  presented 


in  Table  I.  The  average  values  ranged  from  0.3534  to 
0.3862  throughout  the  45  RAPD  loci  and  from  0.2326 
to  0.3093  throughout  the  197  AFLP  loci. 

Population  genetic  structure  -  The  results  of  the 
Wright’s  estimate  (1951)  and  Lynch  and  Milligan’s 
(1994)  F57-and  theta  indices,  using  the  data  from  RAPD 
loci,  ranged  from  0.023  to  0.090;  while  the  AFLP  data 
ranged  from  0.05 1  to  0.224  (Table  II).  These  statistics 
based  on  AFLP  data  were  approximately  double  of  those 
obtained  with  RAPD  data.  Both  markers  showed  the  same 
trend;  that  is,  the  FiT  values  for  population  pairs  using  the 
RAPD  and  AFLP  data,  respectively,  had  higher  values  for 
Choc6-Meta  (RAPD  =  0.055-0.090;  AFLP = 0. 106-0.224) 
and  Cordoba-Meta  (RAPD  =  0.049-0.074;  AFLP  =  0. 1 1 1  - 
0.224)  in  comparison  with  those  for  Choco-C6rdoba 
(RAPD  =  0.008-0.029;  AFLP  =  0.051-0.103). 

Using  the  RAPD  markers,  the  three  gene  flow  esti¬ 
mators  showed  greater  flow  among  the  three  populations 
(Nm  =  2.5-30.4)  than  with  the  AFLP  markers  (Nm  =  0.9- 
4.7)  (  Table  II).  With  both  markers  it  was  observed  that 
there  was  greater  gene  flow  between  Choc6  and  Cordoba. 

Similarly,  the  AMOVA  for  both  sets  of  markers  pro¬ 
duced  higher  <I>ST  values  for  the  analyses  with  AFLP  data 
(Table  III).  The  percent  variation  explained  on  the  basis 
of  the  variance  components  also  showed  appreciable  dif¬ 
ferences  between  the  two  types  of  markers:  8.4  and  23% 
for  RAPD  and  AFLP,  respectively. 

Despite  there  being  high  values  of  intrapopulation 
variation  (92  and  77%  for  RAPD  and  AFLP,  respectively), 
a  significantly  high  interpopulation  variance  was  detected 
for  both  markers  (p  <  0.001).  The  <I>ST  calculated  with 
the  AFLP  data  set  was  similar  to  the  other  F ST  statistics 
calculated  (Table  IV).  On  the  other  hand,  the  estimators 
<I>sr  between  population  pairs  and  the  modified 
coancestry  coefficient  showed  once  again  that  the  popu¬ 
lations  from  Choc6  and  Cdrdoba  had  greater  genetic  dis¬ 
tances  with  respect  to  the  population  from  Meta 

In  confirmation  of  the  foregoing,  the  alleles  frequen¬ 
cies  were  consistently  different  for  a  large  number  of 
loci  from  evaluated  populations.  When  calculating  Nei’s 
(1972)  genetic  distances  between  populations  using  both 
data  sets,  with  a  bootstrapping  of  100  pscudoreplicaies, 
the  consensus  tree  showed  a  greater  genetic  similarity 
between  the  populations  from  Choco  and  Cordoba  than 
between  Mela  and  Choc6  or  Meta  and  Cordoba  (Fig.  3). 

DISCUSSION 

The  results  based  on  the  cluster  analyses,  especially 
for  RAPD  loci,  showed  that  the  Colombian  An.  darlingi 
populations  are  panmitic.  Similarly,  based  on  the  infor¬ 
mation  of  five  RAPD-PCR  primers,  Manguin  ct  al. 
(1999),  demonstrated  evidence  of  geographic  partition¬ 
ing  among  the  analyzed  populations,  however,  the  genetic 
distances  separation  were  very  small  and,  they  concluded, 
that  the  An.  darlingi  populations  are  conspecific. 

The  values  of  genetic  diversity  obtained  in  this  study 
are  nearly  three  times  greater  than  those  observed  by 
Manguin  et  al.  (1999)  (0.063-0.122),  who  used  31 
isozyme  loci  from  seven  populations  from  throughout 
the  range  of  distribution  of  An.  darlingi.  However, 
Mendes  dos  Santos  et  al.  (1999).  working  with  isozymes. 
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TABLKIV 

Distances  <I>ST  among  pairs  of  Anopheles  darlingi  Colombian  populations  with  random  amplified  polymorphic-polymerase  chain 
reaction  (RAPD-PCR)  and  amplified  fragments  length  polymorphism  loci  (AFLP) 


Marker 

Population  pairs 

‘hjj- value  ° 

Modified  co- ancestry 
coefficient  (4>ST) b 

Cordoba  -Mela 

0.1072 

0.1134 

RAPD-PCR 

Choco-Meta 

0.0960 

0.1009 

Cordoba -Choco 

0.0130 

0.0131 

Cordoba-  Meta 

0.2443 

0.2801 

AFLP 

Choco-Meta 

0.2765 

0.3237 

Cordoba-ChtJcd 

0.0860 

0.0899 

a :  distances  <J>yf  among  population  pairs;  b :  distances  tbjy  among  population  pairs,  modified  co-ancestry  coefficient  (-ln|  I  -  «I>ST|=t/2N). 


META 


A 


CHOCO 


1 - 98.0 


CORDOBA 


META 


1 -  100.0 

1 - CORDOBA 

Fig.  3:  genetic  distance  (1-M)  tree  for  three  Colombian  Anopheles  darlingi 
populations  based  on  analysis  of  A:  45  polymorphic  random  amplified 
polymorphic  DNA-polymcrasc  chain  reaction  loci;  B:  197  amplified  frag¬ 
ment  length  polymorphism  loci.  The  branch  number  indicates  the  times 
that  the  topology  was  consistent. 

recorded  higher  values  with  a  greater  observed  and  ex¬ 
pected  intralocus  heterozygosity,  respectively,  in  the 
Manaus  population  (Amazonas)  (Ho  =  0.432  ±  0. 1 1 ;  He 
=  0.375  ±  0.08)  than  for  the  population  from  Cachoeira 
Porteira,  Pani  (Ho  =  0.236  ±  0.09;  He  =  0.290  ±  0.1 1). 
Although  these  values  are  similar  to  those  obtained  by  us 
with  the  aforementioned  markers,  they  are  lower  than  those 
obtained  using  co-dominant  markers  (Yan  et  al.  1999). 

With  data  from  RAPD  loci  of  Aedes  aegypti  from 
Puerto  Rico,  Apostol  et  al.  (1996)  found  an  expected 
heterozygosity  of  0.354,  similar  to  that  found  by  Posso 
el  al.  (2003)  in  An.  nuneztovari  from  Colombia  and  by 
us  in  this  study.  Yan  et  al.  (1999)  discuss  the  variations 
that  can  occur  when  calculating  genetic  diversity  depend¬ 
ing  on  the  marker  type  used.  In  Ae.  aegypti  from  Trinidad 
and  Tobago,  they  found  that  the  heterozygosity  observed 
with  the  RFLPs  was  significantly  higher  (0.47-0.60)  than 
expected  with  the  AFLP  data  (0.39),  values  considered 
to  be  fairly  high. 

The  Fjj  calculations  of  Manguin  et  al.  (1999)  in  the 
seven  An.  darlingi  populations,  obtained  on  the  basis  of 
3 1  isozymatic  loci,  arc  comparable  with  those  for  RAPDs 
(0. 1 02),  but  lower  than  those  for  the  AFLPs  in  this  study. 
Based  on  the  Nm  calculations  between  pairs  of  populations, 
greater  gene  flow  was  observed  between  the  Choco  and 
Cordoba  populations,  as  expected  from  the  cluster  analy¬ 
ses  described  previously;  but  that  was  apparently  overesti¬ 


mated  with  the  RAPD  technique,  given  that  according  to  its 
loci,  the  Nm  values  between  Choc6-Cordoba  were  high, 
especially  based  on  Lynch  and  Milligan’s  F^  (1994). 

The  differences  observed  in  the  Fsr’  estimated  with 
both  types  of  markers  in  this  study,  could  be  explained 
in  the  same  sense  that  Yan  et  al.  (1999)  argue.  They  indi¬ 
cate  that  the  differences  found  among  the  Fsr  estimates 
in  Ac.  aegypti,  when  using  RAPD,  AFLP,  RFLP,  and  iso¬ 
zyme  markers,  could  be  the  result  of  differences  in  the 
mutation  rate  of  these  loci,  which  can  be  examined  in 
Wright’s  (1951)  FST  formula;  Fsr=  l/[l+4(Ve(/n+«)|, 
where  Ne  is  the  effective  size  of  the  population,  in  the 
migration  rate  and  u  the  mutation  rate.  Thus  the  FiT  value 
can  be  seriously  underestimated  if  the  loci  mutation  rates 
are  is  relatively  high.  RAPD  loci  tend  to  have  higher  mu¬ 
tation  rates  than  the  AFLP  and  RFLP  loci;  thus  the  FST 
calculated  from  RAPD  markers  could  be  underestimated. 
On  the  other  hand,  the  isozymatic  markers  mutation  rate 
is  generally  lower  than  for  the  DNA  markers  so  higher 
fS7  estimates  can  be  expected.  According  to  Slatkin  and 
Barton  (1989),  the  theta  statistic  tends  to  overestimate 
the  migration  rate;  but  with  our  RAPD  data  we  observed 
that  although  theta  was  slightly  higher  than  Wright’s  F^, 
it  was,  in  any  case,  threefold  less  than  that  of  Lynch  and 
Milligan  (1994)  for  the  Choc6-Cordoba  comparison.  In 
the  case  of  Lynch  and  Milligan’s  AFLP,  theta  and  F^  the 
data  were  practically  identical;  the  highest  value  was 
obtained  with  Wright’s  FiT 

In  the  AMOVA  (Table  III),  we  observed,  that  the  per¬ 
centage  of  the  variance  components  among  populations, 
was  less  with  RAPD  than  AFLP  markers.  Considering 
that  the  sample  size  used  with  both  markers  types  were 
approximately  the  same,  the  results  differences  should 
stem  basically  from  the  mutation  rate  and  the  number  of 
analyzed  traits  (bands)  (Yan  et  al.  1999).  Based  on  the 
AFLP  loci  analyses,  it  was  observed  that  there  was  an 
indication  of  geographic  subdivision  between  the  west¬ 
ern  and  eastern  populations,  the  same  as  for  the  cluster¬ 
ing  methods.  Between  the  two  western  populations 
(Choco-Cordoba),  there  were  indications  of  a  difference 
comparable  to  random  mating.  For  both  types  of  data, 
there  was  a  high  value  in  the  intrapopulation  variance 
component  (RAPD  =  0.91,  AFLP  =  0.77). 

If  we  consider  that  the  populations  from  Choco  and 
Cordoba  correspond  to  one  ecoregion  and  Meta  to  a  dif¬ 
ferent  one,  the  <I>ST  values  could  correspond  to  a  struc- 
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ture  type  similar  to  that  reported  by  Conn  et  ai.  (1999). 
based  on  mitochondrial  DNA  data  from  samples  of  An. 
darling i  in  Bolivia,  Brazil,  and  Venezuela. 

In  conclusion,  it  is  possible  that  the  eastern  and  west¬ 
ern  populations  of  An.  darlingi  are  genetically  differ¬ 
ent.  Nevertheless,  given  the  delected  differences  with 
the  two  markers  used,  the  search  for  the  population  struc¬ 
ture  should  be  re-analyzed,  using  other  molecular  mark¬ 
ers  such  as  microsatcllites  and  mitochondrial  DNA.  At 
the  same  lime,  parameters  should  be  measured  to  per¬ 
mit  the  analysis  of  its  dispersion,  its  vectorial  capacity 
and  competence,  as  well  as  other  factors  that  may  con¬ 
tribute  to  defining  whether  the  populations  can  be  con¬ 
sidered  cpidemiologically  different. 
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